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ABSTRACT 

The [Nc n] fine-structure emission line at 12.8/im has been detected in several young 
stellar objects (YSO) spectra. This line is thought to be produced by X-ray irradiation 
of the warm protoplanetary disk atmospheres, however the observational correlation 
between [Ne n] luminosities and measured X-ray luminosities shows a large scatter. 
Such spread limits the utility of this line as a probe of the gaseous phase of disks, 
as several authors have suggested pollution by outflows as a probable cause of the 
observed scatter. In this work we explore the possibility that the large variations in the 
observed [Ne n] luminosity may be caused instead by different star-disk parameters. 
In particular we study the effects that the hardness of the irradiating source and 
the structure (flaring) of the disk have on the luminosity and spectral profile of the 
[Ne n] 12.8 fim line. We find that varying these parameter can indeed cause up to 
an order of magnitude variation in the emission luminosities which may explain the 
scatter observed, although our models predict somewhat smaller luminosities than 
those recently reported by other authors who observed the line with the Spitzer Space 
Telescope. Our models also show that the hardness of the spectrum has only a limited 
(undetectable) effect on the line profiles, while changes in the flaring power of the disk 
significantly affect the size of the [Ne n] emission region and, as a consequence, its line 
profile. In particular we suggest that broad line profiles centred on the stellar radial 
velocity may be indicative of flat disks seen at large inclination angles. 
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1 INTRODUCTION 

Young solar mass stars are surrounded by ~ 1O~ 2 M0 of dust 
and gas distributed in a circumstellar disk which is the birth- 
place of planetary systems. In recent years efforts have been 
made to study the inner regions (< 10 AU) of these disks 
in order to constrain gaseous planetary system theories, and 
to understand the physical processes that characterise the 
warm and chemically active regions of protoplanetary disk 
atmospheres. Accurate theoretical modelling of the spectral 
energy distribution (SED) of young stars (see reviews by 
Dullemond et al. 2007 and Natta et al. 2007), coupled with 
observations in the near and mid-infrared carried out with 
the Spitzer Space Telescope have greatly improved our un- 
derstanding of the dust component in the inner disks. On 
the other hand, direct studies of the gaseous component of 
these disks are more challenging, even though the gas ex- 
ceeds the dust by over two orders of magnitude in mass. So 
far molecular line emission, like the CO fundamental and 



overtones (Najita et al. 20071, the rovibrational transition 
of water ([Carr et al.||2008| |Salyk et al.||2008| ), and H 2 pier 



czeg et al. 20021, have been used as diagnostics of regions 



with temperatures around 1000-2000 K. 

The [Ne n] 12.8/^m line (from now on [Ne n]) has been 
suggested as a new probe of gas in the potentially planet 
forming regions of the inner disks (e.g. Glassgold et al. 2007, 
Hollenbach & Gorti, 2009). The ionised region of the disk 
where this line is formed is heated mainly by X-rays from 
the central star (Glassgold et al. 2007, Meijerink et al. 2008, 
Ercolano et al. 2008a, Glassgold, Ercolano & Drake, 2009; 
but see also Gorti & Hollenbach, 2008). Indeed, observa- 
tions from space, with the Spitzer Infrared Spectrometer 
(IRS) JPascucci et al.|[2007l |Lahuis et aLlpOOT] |Espaillat| 



et al.||2007[ ), and from the ground, with MICHELLE and 
TEXES at Gemini North (Herczeg et al., 2007; Najita et 
al., 2009) and VISIR at VLT (Van Boekel et al. 2009), have 
successfully detected the [Ne n] line in tens of young sys- 
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terns. The mean observed line luminosity is comparable to 
the predicted value from disk models, but the scatter of the 
observed values around this mean spans over three orders of 
magnitude (e.g. Giidel et al. 2009). If the [Ne n] line is to 
be used as a useful diagnostic of the gaseous phase of the 
inner disks its origins must be well understood. In particular 
if the line is formed in the disk by irradiation of the gas by 
X-rays from the central star then one naturally expects to 
find a correlation between line luminosity and X-ray lumi- 
nosity for the observed sources. Such correlation is indeed 
observed (Pascucci et al 2007, Giidel et al. 2009), however, 
for objects with the same X-ray luminosity the [Ne n] line 
luminosities span over at least one order of magnitude even 
after removing those sources with known outflows, for which 
it can be argued that a large contribution to the line comes 
from the outflow itself (e.g. T Tau, van Boekel et al. 2009). 
Understanding the origin of this scatter is extremely impor- 
tant if useful information is to be extracted from future ob- 
servations of this line. For example, one should expect that 
disk structural properties, like presence of holes or different 
degree of flaring, not considered in earlier models have an 
effect on the [Ne n] emitting region. Moreover, one relevant 
question is whether the [Ne n] line originates in the bound 
layers of the disk irradiated mainly by X-rays or whether 
it includes a contribution from a photoevaporative outflow 
(e.g. Alexander, 2008). In this paper we will attempt to an- 
swer the latter question by exploring the possibility that 
the scatter observed in [Ne n] line luminosities for a given 
X-ray luminosity may be due to details of the irradiating 
field or to the disk structure (i.e. degree of flaring). We will 
also produce theoretical line profiles and will compare them 
to available high resolution observations (e.g. Herczeg et al. 
2007, Najita et al. 2009). 

We use the models of Ercolano et al. (2009) to calcu- 
late [Ne n] line luminosities and line profiles for irradiating 
spectra of varying hardness. We also show the results from 
additional models that were run to investigate the effects of 
disk flaring on the line luminosities and profiles. The paper 
is organised as follows. The modeling strategy is outlined in 
Section 2; Section 3 contains a review of our results, while 
our final conclusions and discussion are given in Section 4. 



2 MODELLING STRATEGY AND METHODS 

We investigate the effects of two parameters on the [Ne n] 
emission region of protoplanetary disks: (i) hardness of the 
irradiating spectrum; (ii) flaring of the disk. 



2.1 Hardness of the irradiating spectrum 

We calculated fine structure [Ne n] line luminosities using 
the temperature and ionisation structure obtained by Er- 
colano et al. (2009) with the MOCASSIN code for EUV+X- 
ray irradiated protoplanetary disks. We refer to Ercolano et 
al. (2008a, 2009) for details of the disk models, and to Er- 
colano et al. (2003, 2005, 2008b) for details of the code used. 
The atomic database included opacity data from Verner 
et al. (1993) and Verner & Yakovlev (1995), energy lev- 
els, collision strengths and transition probabilities from Ver- 
sion 5.2 of the CHIANTI database (Landi et al. 2006, and 
references therein) and the improved hydrogen and helium 
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Figure 1. Isodensity curves of the two input disk models used 
in this work to explore the impact that different degree of flaring 
has on [Ne n] line luminosity. The black solid lines show the 'flat' 
model (/3 = 1.0) and dashed lines show the 'flared' model (/3 = 
1.25). The plotted contours refer to increasing densities from 10 5 , 
at higher heights, to 10 13 g cm -3 deeper into the disk with steps 
of 10 2 cm" 3 . 



free-bound continuous emission data of Ercolano & Storey 
(2006). In particular the Ne + atomic data include collisional 
strengths of Saraph & Tully (1994) and the transition prob- 
abilities of Blackford & Hibbert (1994). We briefly sum- 
marise here the model input parameters: we used a stellar 
mass of 0.7 Mq, radius 2.5 Kq and effective temperature 
of 4000 K and an initial disk mass of 0.027 Mq with an 
outer radius of ~500ATj(3 The MOCASSIN models are cou- 
pled to a hydrostatic equilibrium routine which ensures that 
the effects of irradiation on the disk structure are taken 
into account self-consistently. Line luminosities computed 
without such routine result to be underestimated at worst 
by a factor of ~ 2. The models were irradiated with syn- 
thetic spectra extending from the EUV to the X-ray spec- 
tral region (13.6 eV-10 keV) where the spectra were calcu- 
lated using the multi-temperature plasma prescription de- 
scribed by Ercolano et al. (2009) and were computed us- 



ing the PINTofALE IDL software suite ( Kashyap & Drake 
|2000[ ), adopting the solar chemical composition (Grevesse & 
Sauval, 1998), the atomic data from the CHIANTI compila- 
tion (Landi et al. 2006) and ion populations from Mazzotta 
|et al.| ( |1998| ). The transmittance of the synthetic spectra 
through neutral screens of varying thickness (as expected to 
be found in the circumstellar environment of pre-main se- 
quence stars, e.g. Giidel et al., 2007, 2008) was calculated 
assuming neutral hydrogen column densities in the range 
Nh = 10 18 — 10 22 cm -2 , with increasing columns result- 
ing in increasingly harder spectra, as the soft radiation is 
absorbed out. We present [Ne n] luminosities and line pro- 
files for the Ercolano et al. (2009) models labeled in their 
work as FS0H2Lxl, FS18H2Lxl, FS19H2Lxl, FS20H2Lxl, 
FS21H2Lxl and FS22H2Lxl, which consist of models ir- 
radiated by spectra screened by columns of 0., 10 1S , 10 19 , 
10 20 , 10 21 , and 10 22 cm , and an X-ray luminosity after the 
screen of Lx(0.1-10 keV) = 2xl0 30 erg/s ("base luminosity" 



1 Note that only the inner 50 AU of the disk are considered here 
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hereafter; please refer to Table 1 in Ercolano et al. (2009) for 
a full description of the model parameters). Keeping Lx (0.1- 
10 keV) constant after the screen by normalization of the ab- 
sorbed spectrum was done here in order to isolate the effects 
of the incoming radiation field hardness alone. Our prescrip- 
tion results in the unscreened models to have a total ionis- 
ing luminosity L to t(13.6 eV-10 keV) of approximately twice 
the X-ray luminosity and the moderately screened models 
(Nh = 10 18 cm -2 ) to have a total ionising luminosity 1.2 
times the base luminosity. We ran four additional models 
with screens of and 10 21 cm -2 and X-ray luminosities an 
order of magnitude higher and lower than the base luminos- 
ity in order to ensure that any scatter produced by varying 
screen thickness at fixed Lx = 2xW 30 erg/ s would also be 
reproduced for higher or lower X-ray luminosities. 



2.2 Disk Flaring 

In order to investigate the effect of disk flaring on [Ne ll] 
line luminosities, we used a simple analytical description of 
the disk structure, as given by Robitaille et al (2006). The 
density distribution is described by a fixed power-law decline 
of the surface density in the radial direction coupled with a 
gaussian relationship along the z-direction: 



p(R,z) 




exp 



2h 2 (R) 



(1) 



where po is a scale factor fixed by the disk total mass and 
h(R) is the disk pressure scaleheight parametrized by a 
power- law: : 



h{R) 



(2) 



where Ri„ and hi„ are, respectively, the dust destruction 
radius and the hydrostatic equilibrium scaleheight at this 
radius. Different values of /3 produce disk models with dif- 
ferent degree of flaring, with higher (3 producing higher flar- 
ing. We adopted /? = 1.0 and f3 = 1.25 to build models for 
a 'flat' and 'flared' disk, respectively, and we calculated the 
corresponding [Ne n] luminosities and line profiles. Such val- 
ues are also used as limits for the (3 parameter in the set of 



SEDs of Robitaille et al. ( 2006 1, typically used in the litera- 



ture to fit real observations. For either density structure we 
also chose Ri„ corresponding to a dust sublimation tempera- 
ture of 1600 K, i.e. Ri„ ~0.08 AU with the previous central 
star parameters, z factor = 1 from model comparison with 
the disk structure of Ercolano et al. ( 2009 1 and a = (3 + 1 



( Whitney et al.|2003a|b " ;o preserve the surface density dis- 
tribution power law scaling E(R) ~ R~ 1,5 of the minimum 
solar nebula model of Hayashi (1981). The other stellar and 
disk parameters were chosen to be the same as those used 
by the Ercolano et al. (2009) models described above. 

Figure [T] shows the isodensity curves of the two disk 
models with different flaring powers. The flat model has a 
scaleheight reduced typically by a factor of about 3-4 be- 
tween cylindrical radii from 1 to 10 AU when compared to 
the flared model. In order to isolate the effects of disk flaring 
only we keep the density structures fixed and do not impose 
hydrostatic equilibrium on these models. 



2.3 Line profiles 

We computed line profiles from all our models assuming 
Keplerian rotation at several disk inclinations. The profiles 
were computed by summing up the contributions from each 
cell which produces a Doppler broadened profile given by 



I(R,z) 



/2irvth{R, z) 



exp 



{v-vi os {R,e)f 

2v th (R, z) 2 



(3) 



where I(R, z) is the power emitted by each cell and Vth is 
the local thermal velocity of the Ne + ions equal to 



v t h(R, z) 



3k B T(R,z) 



(4) 



and vios {R, 0) is the cell velocity projected along the line 
of sight. Given the temperature structure of the disk, the 
width of the thermal broadening in each cell is ^6 kms" 1 . 
We assumed that the midplane is optically opaque to the 
[Ne ll] line radiation so only contributions from one side of 
the disk are taken into account when computing the profile. 
In fact, adopting the dust cross section from Weingartner 
and Draine 2001 of Rv = 5.5 extinction law, generally used 
to reproduce the interstellar extinction in molecular clouds, 
0-12. s„ m = 1.6 xl0~~ 23 cm 2 / H and tris.s^m = 1.4 x 10~~ 23 
cm 2 / H, we find that optical depths of the order of 1 are 
reached for hydrogen column density of about 6-7 x 10 22 
cm -2 . Such column density are reached, for all inclination 
5j 60, at heights of ^ 0.1 AU at a radial distance of 1 AU 
and ~1 AU at a distance of 10 AU, deeper in the disk than 
the Ne emitting region shown in Figure [3] In the case of 
almost edge-on systems the column density along the line 
of sight could be large enough to attenuate the observed 
line intensities, this is not taken into account in this work, 
but should be in future comparisons with real observations. 
Finally, we ignore scattered light contribution to the profile, 
which is negligible in the Weingartner & Draine (2001) dust 
model that has an albedo of order 0.01 at 12.8 pm, but could 
be significant for different dust models. The resolution of 
the line profiles we computed is A/AA ~ 150000, since we 
summed all together the emission from cells whose vi os falls 
in bins wide 2 kms -1 . Moreover, we have convolved the line 
profiles with a gaussian function with FWHM of 10 kms" 1 
to degrade them to a resolution A/AA ~ 30000, comparable 
with the high-resolution spectroscopic observations. 



3 RESULTS 

In this section we present the line luminosities and profiles 
obtained from the models described above. Monte Carlo er- 
rors on the given line luminosities were obtained by com- 
paring multiple runs of the same model and are about 5% 
for the lines in question. As in Sect. 2. 3 we will consider the 
[Ne ll] line luminosity response to (i) changes in the shape of 
the irradiating spectrum and (ii) changes in the disk density 
distribution (flaring). 



3.1 Hardness of the irradiating spectrum 

Table [j] shows the values of [Ne n] 12.8 /im and [Ne in] 15.5 
/im line luminosities for models irradiated by the multi- 
temperature thermal spectrum of Ercolano et al. (2009) 
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Figure 2. L([Ne ll])/Ltot as a function of the screening column 
of the irradiating spectrum. The dashed line indicates the value 
for the unscreened spectrum. 
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Figure 3. Two-dimensional distribution of the [Ne II] emissiv- 
ity in arbitrary units for a hydrostatic equilibrium disk model 
irradiated by an unscreened spectrum. 



screened by columns of circumstellar material of thickness 
0., 10 18 , 10 19 , 10 20 , 10 21 , 10 22 cm" 2 (see their Figure 2). 

Figure [5] shows L([Ne n])/L to t as a function of the 
screening column. The dashed line indicates the value pre- 
dicted in the case of unabsorbed source. From this figure and 
from the values in Table [I] it appears that EUV radiation 
is less efficient at producing [Ne n] than the (soft) X-rays. 
In fact, even if the unscreened models have a total ionising 
luminosity (X-Rays+EUV) which is double than the X-ray 
only luminosity, they only produce a moderate increase in 
L([Ne n]). This can be easily understood when one considers 
that [Ne n] is most efficiently produced in the low ioniza- 
tion layer (just above the molecular layer, see Glassgold et 
al., 2007; Glassgold, Ercolano & Drake, 2009), where a fast 
charge exchange with the abundant neutral H atoms quickly 
transforms most doubly ionised Ne into singly ionised. The 
charge exchange coefficient from singly ionised to neutral 
Ne is instead extremely slow, conspiring to produce high 



abundance of Ne + in this warm layer. Instead, EUV radi- 
ation (> 21.6eV and < lOOeV) can only produce Ne + by 
standard ionisation of neutral Ne and is only efficient in the 
lower density upper atmosphere where an H II region like 
layer is formed. However, some of the Ne + produced in this 
region will be ionised to Ne 2+ and the lower neutral H frac- 
tions in this region will result in a much less efficient charge 
exchange process. A spatial map of the [Ne n] emission for 
the unscreened model is shown in Figure [3] 

Furthermore, the models shown in Table [I] also suggest 
that the hard X-ray region is inefficient at the production of 
[Ne n] compared with the soft X-ray region. Indeed models 
irradiated by hard spectra with Nh ^ 10 21 cm 2 show a 
steep decline in the [Ne n] luminosity, even after forcing 
the X-ray luminosity (0.1-10 keV) to stay constant. This 
effect is dominated by the fact that models illuminated by a 
harder X-ray spectrum result to be colder (see discussion in 
Ercolano et al. 2009), the hard X-ray penetrate a region that 
is denser and cooler and therefore the resulting [Ne n] flux 
is reduced. Cooler models also result in flatter disks which 
subtend a smaller solid angle to the source of radiation. 

The hardness of the spectrum has a noticeable effect 
on the [Ne in] 15.5 /im/[Ne n] 12.8 /im ratio. This ratio is 
higher for softer spectra due to the fact that doubly ionised 
Ne in these models can be produced in the upper layers of 
the disk where the neutral hydrogen fraction is lower and 
therefore charge exchange less efficient. Observations of the 
[Ne ill] 15.5/im (from now on [Ne ill]) line are still rather 
sparse although so far low values of the [Ne m]/[Ne n] ratio 
have been inferred (e.g. 0.06, Lahuis et al., 2007), indicating 
that according to our models, at least for the few objects 
in question, a soft (EUV-rich) irradiating spectrum is not 
expected. 

Our results show that changes in the irradiating spec- 
trum can give an order of magnitude scatter (an extra factor 
of ~2 is introduced when one considers that an unscreened 
source may contribute comparable ionising luminosity in the 
EUV band that would not be accounted for in the typical 
L([Ne n]) vs Lx relation). 

Figure [5] (left panel) shows the surface-integrated emis- 
sion for some of the models discussed in this Section. We 
find that in all cases half of the total [Ne n] emission comes 
from the inner regions of the disk, < 3-4 AU, so that the dif- 
ferent spectra do not significantly affect the emission region 
for this line. Table [3] lists the half-luminosity and 90% lumi- 
nosity radii for models of various screens. As the extent of 
the emission region does not significantly vary with spectral 
hardness, the line profiles obtained are also roughly invari- 
ant. As an example, in Figure [6] we show the [Ne n] profile 
for the unscreened model seen almost face on (10° inclina- 
tion) and almost edge on (80° inclination), which have Full 
Width Half Maximum (FWHM) of ~6-8 and ~25-31 kms' 1 
respectively. 



3.2 Disk Flaring 

Table [5] also shows the results for the L([Ne 11]) for models 
with flared (/3 = 1.25) and flat (/? = 1.0) disks. A casual in- 
spection of this table reveals roughly an order of magnitude 
increase in the relative luminosity of the [Ne 11] line when 
the flaring is changed from j3 = 1.0 to /3 = 1.25. This is 
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Table 1. [Ne ll] and [Ne III] luminosities of disk irradiated by 
different high-energetic spectrum of Lx(0.1 — lOkeV) = 2xl0 30 
erg s _1 . The disk model adopted is a d'Alessio disk model, see 
also Ercolano et al. 2009 for details. 



Table 3. Emissivity distribution statistics. See text for details. 



MODEL L tot (13.6 eV-10 keV) 

10 30 erg s" 1 



[Ne n] [Ne in] 

erg s —1 ) 



(10 28 



Unscreened 
Log N h = 18 
Log N h = 19 
Log N h = 20 
Log N h = 21 
Log N h = 22 



1.0 
2.3 
2.04 
2.0 
2.0 
2.0 



2.3 
1.9 
1.7 

1.5 
9.1x10" 
1.4x10" 



4.9x10 
4.3x10 
3.2x10 
2.2x10 
1.5x10 
2.5x10 



Model 


Half-L radius 


90% L radius 




(AU) 


(AU) 


Unscreened 


2.9 


10.1 


Log N h = 20 


3.3 


13.0 


- Log N h = 21 


3.5 


12.8 


Log N h = 22 


2.9 


10.1 


Flared (3 = 1.25 


12.9 


34.3 


Flat (3 = 1.0 


1.5 


6.1 



Table 2. [Ne ll] line luminosities computed at a range of Lx ■ The 
corresponding values of Ltot are also given. 



MODEL [Ne n] 

(10 28 erg s" 1 ) 



L x [210 30 ergs -1 ] 


0.1 




1.0 


10. 


Unscreened (Ltot = 2Lx) 


2.9x10" 


-l 


2.3 


22.0 


Log(N H )=21 (Ltot = Lx) 


3.0x10" 


-2 


9.1X10" 1 


5.8 


Flared (L to t = 2L X ) 


1.0x10" 


-1 


1.6 


17.3 


Flat (Ltot = 2L X ) 


1.7x10" 


-2 


1.3X10" 1 


7.0x10 



Table 4. FWHM of the line profiles simulated for disks in Keple- 
rian rotation viewed almost face on (10° inclination) and almost 
edge on (80° inclination) at two different spectral resolutions. 



Inclination 10 degree Inclination 80 degree 
FWHM (kms- 1 ) FWHM (kms" 1 ) 

R = 150000 30000 150000 30000 



Unscreened 8.2 11.1 31.3 32.6 

Log N H = 21 6.2 10.0 25.8 28.1 

Flared (3 = 1.25 6.9 10.2 16.8 17.6 

Flat (3 = 1.0 9.5 12.1 40.1 40.8 



easily understood when considering that a more flared disk 
would subtend a larger solid angle to the irradiation source. 

Disk flaring has also a dramatic effect on the size of 
the [Ne n] emitting region. Figure [4] shows two-dimensional 
plots of the [Ne n] emissivity for the flared and flat disk 
models, while the right panel of Figure [5] shows the surface- 
integrated emission for the line for the two disk models. As 
summarised in Table [3] the emission region in the flared 
disk model extends out to a radius of ~30 AU, within which 
90% of the total [Ne ll] luminosity is emitted. This radius is 
reduced to ~6 AU for the flat disk model. Naturally, this has 
consequences on the line profiles for these two cases, which 
are shown in Figure [7] for two disk inclinations. The FWHM 
for the [Ne ll] line produced by a flat disk seen almost edge 
on (inclination 80°) is roughly 40 km s _1 , while it is only 
17 km s _1 for a flared disk where the emission region extends 
to larger cylindrical radii. The difference is of course smaller 
if the disks are seen almost face on (inclination 10°), thus 
dominated by the contribution of the thermal broadening, 
where a FWHM of ~10 km s"" 1 is expected for a flat disk 
compared to a FWHM of ~7 km s _1 for a flared disk. These 
statistics are summarised in Table [4] 



4 DISCUSSION 

As the sample of observed L([Ne ll]) luminosities from pro- 
toplanetary disks continues to grow, together with the sam- 
ple of pre-main sequence stars with measured X-ray lumi- 
nosities, a trend between the two quantities becomes more 
and more evident (Gudel et al. 2009). While this lends evi- 
dence to the [Ne ll] line being a product of X-ray irradiation 



of disks as predicted by Glassgold et al. (2007), the large 
scatter associated with this relation makes the situation less 
clear-cut. For a handful of objects with known powerful out- 
flows it is clear that this line cannot be used as a diagnostics 
of the gaseous disk phase, as the emissivities may be domi- 
nated by gas in the outflows. However, even after removing 
objects with known outflows from our sample we are left 



1.5 




Figure 6. [Ne ll] line profiles of a Keplerian disk in hydrostatic 
equilibrium irradiated by an unscreened EUV+X-ray source. The 
black lines are for low inclination (i = 10°) and the red lines for 
high inclination (i = 80°) disks. The solid lines have a resolution 
of R = 150000 and the dashed are the profiles degraded to a 
resolution of 30000. For clarity the lines are normalized to the 
peak value of the degraded case. 
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Figure 4. Two-dimensional distribution of the [Ne II] emissivity in arbitrary units. The left panel refers to the disk model with /3 = 1.25 
(flared disk), instead the right panel show the disk model with /3=1.0 (flat disk). 
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Figure 5. [Ne II] cumulative surface emission in erg s _1 (see text for details). Left: Hydrostatic equilibrium models irradiated by spectra 
screened by various neutral hydrogen columns. Right: Flat and flared disk models computed assuming Robitaille et al. (2006) density 
distribution. 




Figure 7. As in Figure|6] but for the flared disk model (0 = 1.25) (left) and fiat disk model (f3 = 1.00) (right) irradiated by an unscreened 
EUV+X-ray spectrum. 
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with approximately one order of magnitude variations in 
L([Ne ll]) at any given Lx- 

In this paper we have investigated the origin of the large 
scatter observed in the [Ne ll] - Lx relation. We find that 
variations in the irradiating spectral shape and disk struc- 
ture (flaring) are sufficient to explain the typically observed 
scatter of approximately one order of magnitude. Figure [8] 
shows the results from our models (coloured symbols) over- 
plotted on a shaded region representing roughly the observa- 
tions presented by Giidel et al. (2009), excluding the strong 
outflow objects (blue points in their plot). The definition of 
such region is affected by the limited sample statistic that 
does not allow a sharper outline, but it shows where the bulk 
of data falls. Moreover, we note here that we use luminosities 
after the screen while Giidel et al. use intrinsic luminosities 
derived from the observations. We have performed absorp- 
tion calculations for a subset of the observations from which 
we found the following: absorbing the intrinsic X-ray spec- 
tral models by N H = 3 x I0 20 cm -2 - 4.6 x 1 22 cm~ 2 (a 
range occupied by the observations) reduces the 0.1-10 keV 
luminosity by factors of approximately 1.4-5.5 (also depend- 
ing on the intrinsic thermal model from the individual fits). 
On average, thus, the unabsorbed Lx from the observations 
should be shifted by a factor of 3 (0.5 dex) to the left to ob- 
tain average absorbed values, with individual shifts ranging 
from 0.15 dex to 0.75 dex. These corrections do not change 
the plot qualitatively. We, however, refrain from using the 
individual absorbed (post-screen) luminosities derived from 
the observations here because it is not clear a priori that 
the absorbed luminosity seen from Earth is close to the ab- 
sorbed luminosity seen by the disk. The latter luminosity 
is principally unknown. We also mention that Giidel et al. 
use the 0.3-10 keV range for their luminosities, while we use 
0.1-10 keV. For the post-screen luminosities for stars with 
Nh > 10 20 cm -2 , however, the difference is minor. 

Our model predictions seem to agree well with the ob- 
servations in terms of reproducing the scatter, and at in- 
termediate and high X-ray luminosities the absolute val- 
ues of L([Ne ll]) also roughly agree with the observations. 
Models with low X-ray luminosities, however, fall short of 
the observed values, and while detection limits can par- 
tially affect the observations, the failure of our models to 
produce L([Ne n]) > I0 28 erg/ sec for X-ray luminosities of 
2xl0 29 ergs~ 1 must have another origin. Nevertheless there 
is an overlap between our prediction for the unabsorbed 
source and/or the flared disk and the shaded region. Keep 
in mind, however, that the correction for absorption men- 
tioned above aggravates the discrepancy between the model 
predictions and the observations (the latter being moved 
somewhat to the left in Fig. 8), but do not have effects on 
the conclusions on the scatter of [Ne ll] luminosities. 

It would be of little use here to speculate what other 
observational effects may be coming into play, however one 
thing that is worth noticing is that the measurements plot- 
ted in the Figure all come from the Spitzer Space Telescope, 
while recent observations from the ground at higher resolu- 
tion suggest often lower luminosities than the Spitzer data 
( jHerczeg et al.||2007| [Najita et aL|2009| >. This has been in- 
terpreted as possible pollution from molecular emission in 
the Spitzer band, while more recent works (Najita et al., 
2009) seem to lean toward extended emission, from unde- 
tected outflows, which would contribute to the flux in the 



Spitzer aperture, but would be excluded by the narrow slits 
used for measurements from the ground. 

Our results thus support a disk origin of the observed 
[Ne ll] emission as a consequence of X-ray irradiation, at 
least for systems with moderate [Ne ll] luminosities and ab- 
sence of jets. On the other hand, the additional parameters 
(X-ray spectral hardness, disk flaring) now found to influ- 
ence the [Ne ll] luminosity add complexity to the interpreta- 
tion, and the usefulness of the [Ne ll] line as a disk diagnostic 
depends on our ability to confine these other parameters as 
well. Clearly, [Ne ll] disk diagnostics needs to avoid stars 
with known jets and outflows. Under the assumption of full 
dust and gas mixing, some information about disk flaring 
can be obtained by modeling the infrared spectral energy 
distribution ( Chiang et al.|[2001 Pascucci et al.|2 003), pro- 
viding some important constraints on the range of modeled 
[Ne ll] luminosities. Constraining the shape of the ionizing 
spectrum is perhaps the most difficult task. The intrinsic 
spectrum can usually be modeled sufficiently well based on 
intermediate- or high-resolution X-ray spectroscopy, but the 
absorbing gas column density toward the disk may be differ- 
ent from that toward the observer. However, gas absorption 
matters the least in systems with weak or absent accretion, 
the possible presence of inner holes, and therefore they likely 
absence of strong winds emanating from the innermost disk 
regions. The [Ne n] disk diagnostic may thus be optimized 
for the study of transitional disks, most of which are also 
not known to drive outflows or jets. Several of these have 
been well studied, including new ground-based observations 
of spectrally resolved [Ne ll] lines (Najita et al. 2009). In 
this context, the [Ne ll] line may develop its full diagnostic 
power to study disk ionization and heating by ionizing ra- 
diation from the central star, processes that are pivotal for 



disk dispersal through photoevaporation (Alexander et al 



2006 Ercolano et al. 2008 2009 Owen et al. 2009) in the 



first place. 

High resolution and high signal-to-noise line profile can 
help assessing the origin of the line, with profiles centred at 
the stellar radial velocity being synonymous with a disk ori- 
gin. This is true for all inclinations except those very close 
to edge on, where outflows would also produce profiles cen- 
tered on the stellar radial velocity. Blue shifted profiles are 
expected to be observed for system with outflows observed 
at non edge on inclinations (e.g. Alexander, 2008), and have 
been recently oberved by Pascucci et al. (2009) in the spec- 
tra of TW Hya, Sz 73, T Cha and CS Cha. Other examples 
of high resolution spectra where the line has been detected 
include TW Hya (Herczeg et al. 2007; A/ A A ~ 30000), GM 
Aur and AA Tau (Najita et al., 2009; A/AA ~ 80000). Here 
the observed lines are consistent with being centred at the 
stellar radial velocity, suggesting a circumstellar disk origin, 
although the poor signal-to-noise of some of these detections 
makes it hard to say with certainty. The lines are broad with 
a FWHM of respectively ~ 21 (TW Hya), 70 (AA Tau) and 
14 (GM Aur) km s^ 1 . The [Ne ll] FWHM of GM Aur can 
be produced by a disk of normal flaring power, however the 
large FWHM of AA Tau requires the emission region to be 
dominated by very small radii and therefore implies a small 
degree of flaring. We finally conclude that high resolution 
observations of the [Ne n] line in YSOs, able to resolve its 
profile, are needed if this line is to be used to extract useful 
information on disk structure (e.g. flaring) and evolutionary 
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Figure 8. Comparison between the prevision of our models, 
where a relationship between the L([Ne n]) and hx is expected, 
and the observed system with undetected outflows (from Giidel 
2008). 



stage by (e.g.) the detection of outflows and photoevapora- 
tive winds by comparison with line profile models like those 
of Alexander (2008) or those presented in this work. 
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